Light self-trapping in a large cloud of cold atoms 
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We show that, for a near-resonant propagating beam, a large cloud of cold 87 Rb atoms acts as a saturable 
Kerr medium and produces self-trapping of light. By side fluorescence imaging we monitor the transverse size 
of the beam and, depending on the sign of the laser detuning with respect to the atomic transition, wc observe 
self-focusing or -defocusing, with the waist remaining stationary for an appropriate choice of parameters. We 
analyze our observations by using numerical simulations based on a simple 2-level atom model. © 2013 
Optical Society of America 
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In 1974, Bjorkholm and Ashkin observed the first opti- 
cal spatial soliton in a bulk nonlinear medium, namely 
a hot sodium vapor [1]. Indeed, it was previously ar- 
gued [2] that the diffraction associated to a finite trans- 
verse size could be balanced by self-focusing e.g. in a 
Kerr medium, where the intensity profile of the laser 
beam induces a positive lens. This self-trapped beam 
can propagate in a quite stable fashion if the nonlinear- 
ity is saturable. Spatial solitons have, since then, been 
studied in a variety of nonlinear media [3], including, 
for instance, glass [4], photorefractive crystals [5], liq- 
uid crystals [6] and light- valves [7]. In the last decade, 
laser-cooled atomic vapors have become increasingly em- 
ployed to study both linear and nonlinear light propaga- 
tion. Compared to other media, cold atoms offer several 
advantages such as e.g. the possibility to manipulate the 
properties of the medium, the large nonlinearities at low 
light levels, and the simple direct modeling from the mi- 
croscopic expression of the atomic susceptibility. A few 
preliminary studies [8-10] have addressed the issue of 
self-focusing and -defocusing in cold atomic samples, in 
the "thin-lens" regime where modifications of the beam 
profile inside the atomic medium are negligible. 

In this Letter, we present the first experimental evi- 
dence of light self-trapping inside a cloud of laser-cooled 
atoms. This is achieved by sending a near resonant beam 
to propagate inside a large cloud of cold 87 Rb atoms. 
The solitonic behavior is detected by monitoring the size 
of the beam through fluorescence side-imaging. After a 
brief description of the setup, we present our observa- 
tions on the self-focusing and -defocusing regime and we 
show that there exists a regime of parameters for which 
the beam size remains stationary. The results are sup- 
ported and discussed with the help of numerical simula- 
tions. 

The experimental setup is sketched in Fig. 1. Observ- 
ing self-trapping requires the length of the nonlinear 
medium to be quite larger than the Raylcigh length of 
the incoming laser beam. Moreover, the spatial atomic 
density should be large enough to yield a sizable index of 
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Fig. 1. (Color online) Experimental setup. We focus a 
laser beam (waist 20 /im) at the front side of the cold 
atomic cloud, and monitor its transverse size as it prop- 
agates using side fluorescence imaging. The far-field dis- 
tribution of the transmitted beam is also recorded. 



refraction. Fulfilling both criteria with cold atoms repre- 
sents quite a challenge, requiring trapping and cooling a 
large number of atoms. We prepare our sample of cold 
atoms in a "giant" vapor-loaded MOT using six trap- 
ping beams of waist 3 cm with a total power of 600 mW, 
and containing roughly 10 11 atoms at a temperature of 
« 100/xK. Because of light-induced inter-atomic repul- 
sion, the cloud size is large with a quasi-Gaussian den- 
sity profile of peak value po — 8 x 10 10 cm -3 . Along the 
beam propagation axis z, the FWHM size of the cloud 
is 11 mm, much larger than the probe's Raylcigh length 
Zr = 1.6 mm (beam waist Wi n = 20 yum at the front 
entrance of the cloud). Once the MOT is loaded (within 
Is), we switch off the trapping lasers and magnetic field 
and shine the probe beam at the atoms for a 2 (is dura- 
tion. This short time is chosen to avoid nonlinear effects 
linked to optical pumping or mechanical effects [10]. We 
checked that increasing the pulse duration up to 10 /xs 
does not affect our observations. In these conditions an 
atom scatters at most 38 photons during the interaction 
with the laser beam. We record from the side the fluo- 
rescence profile induced by the beam, using an EMCCD 
(Andor iXon+). The image is accumulated for up to 100 
cycles. The far-field distribution of the transmitted beam 
is recorded using another CCD camera. 

Fig. 2 shows an evidence of self-trapping of the laser 
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beam inside the atomic cloud. We plot the transverse 
FWHM size Ay (as obtained from a Gaussian fit) of 
the fluorescence induced by the probe, as a function of 
z. We compare the propagation for opposite detunings 
from the atomic resonance 6 = lji — Lu at , where uii and 
oj a t are the laser and atomic angular frequencies, respec- 
tively. The solid and open circles correspond to 6 = 5T 
and — 5r respectively, where T = 2tt x 6.06 MHz is the 
natural lincwidth. The laser power is 100 fiW. The meas- 
ured transmissions are 26% for 5 > and 18% for <5 < 0. 
The bold curves are obtained from numerical simulations 
(see below) whereas the curve at bottom is the measured 
atomic density profile. The observed behavior clearly de- 
parts from linear propagation and shows a strong de- 
pendence on the sign of the detuning. For 6 > 0, after 
an initial growth the beam size remains constant, close 
to the expected behavior for a self-trapped beam when 
self-focusing exactly cancels diffraction. This behavior 
is emphazised in the inset, where we compare situations 
where self-focusing overcomes diffraction (8 = 3T, empty 
circles), balances it (5 = 5r, solid circles) or is dominated 
by it (5 = 7r, empty squares). 

From Fig. 2 we see that self-focusing takes place on 
typical distances of the order of the cloud size, hence 
the thin-lens description [8-10] cannot be used here. A 
quantitative analysis requires full numerical simulations 
of the nonlinear propagation inside the cloud. It relies on 
a 2-level model, valid in the absence of optical pumping 
(Zeeman or hyperfinc) and mechanical effects [10], which 
is justified for our short pulse duration. In the absence of 
Doppler broadening, the refractive index n and photon 
scattering rate per atom V are given by 

tr n i 3A3 s JL m 

n ^ = 1 -^i + m rf + 9 i/i sat (1) 
// r gi/isat v s 

r(/ ' } = 2 1 + A(5/Yf + gl/l sat = 2 TT^ (2) 

where s(I,S) = g(I/I sat )/[l + 4(S/T) 2 ] is the satura- 
tion. Here, g — 7/15 is the degeneracy factor account- 
ing for the presence of several equally-populated Zee- 
man sub-levels, p the atomic density, I the laser inten- 
sity and I S at = 16 W/m 2 the saturation intensity. The 
dispersive behavior of the refractive index (Eq. 1) yields 
self- focusing and -defocusing for 6 > and < 0, respec- 
tively. The photon scattering rate (Eq. 2) is responsible 
for the attenuation of the laser beam as it propagates 
through the sample. The measured fluorescence signal 
is S(y,z) oc p(z) fT'(x,y,z)dx, where we assume that 
the transverse atomic density profile is uniform on the 
scale of the beam size. Thus, the fluorescence signal does 
not directly reflect the "beam size" I(y, z), except in the 
linear case s -C 1. In addition, the scattered photons 
generated in the laser beam have to propagate through 
the atomic cloud to reach the CCD. This propagation 
is linear outside the beam, but diffusive if the optical 
thickness of the cloud along x is > 1. Assuming scattered 
photons at the laser frequency, this condition is fulfilled 



for |<5| > 3r. Indeed, we typically observe a diffusive 
pedestal on our fluorescence beam profiles for \S\ < 4T. 
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Fig. 2. (Color online) Observation of self-trapping. We 
plot the fluorescence transverse size Ay versus position z 
in the atomic cloud. We compare symmetric detunings: 
S = 5r (solid circles) and S = — 5T (open circles). The 
beam waist is 20 fim, positioned at z ~ —11.5 mm, and 
the laser power is 100 /LtW. The bold lines are numer- 
ical simulations (see text). The curve at bottom is the 
measured atomic density profile. The inset compares the 
measured data for 6 = 3T (open circles), 5 = 5T (filled 
circles) and 5 — 7T (open squares). 

A qualitative picture for self-trapping can be provided 
assuming a "step-like" transverse profile for the refrac- 
tive index. This assumption is reasonable for a uniform 
circular illumination [2] , or for a Gaussian beam at high 
saturation s 3> 1 [10]. If the refractive index at beam 
center is larger than that in the wings, which is the case 
for 5 > (Eq.l), a critical angle 9 C for total internal re- 
flexion can be defined. Thus, self-trapping can occur if 
the incident beam divergence 6d = X/nu)i n satisfies 



Srf< r^f^iTlrF' (3) 

This inequality is independent of / since we assumed 
s 3> 1. For Win = 20 /urn and 8 = 5T, self-trapping 
occurs for p > 8 x 10 10 cm -3 . In addition to assuming 
a step-like profile (s 3> 1), this simple approach neglects 
two important features of our system. First, the density 
profile (and thus the index step) is not uniform along the 
propagation direction. This implies that the condition 
above can be satisfied only in a limited spatial range 
around the center of the cloud, which explains the initial 
increase of the beam size for S > in Fig. 2. Second, with 
our parameters the attenuation of the beam is rather 
large due to scattering. This attenuation limits the range 
of z where the condition s>l can be sustained. 

A numerical simulation is performed, including the 
non-uniform density profile of the medium along z and 
the beam attenuation due to photon scattering. We em- 
ployed a split-step method, computed the intensity dis- 
tribution I(x, y, z) and, then, the fluorescence signal and 
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the far-field distribution of the output beam. The bold 
curves in Fig. 2 show the behaviors obtained with the 
measured experimental parameters, except for the peak 
atomic density (1.1 x 10 11 cm" 3 here) which was in- 
creased by a factor 1.3 to compensate for a systematic 
discrepancy between experiments and simulations. Such 
an uncertainty of 30% on the peak density is reasonable, 
since this quantity is difficult to determine accurately 
for large MOTs because of high optical thicknesses and 
non-Gaussian profiles. The agreement is then quite good 
both qualitatively and quantitatively. 
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Fig. 3. (Color online) Nonlinear propagation versus laser 
power. The divergence (see text) is plotted vs the laser 
power, for S = AT (circles) and 5 = —4Y (squares). The 
filled and empty symbols correspond to experiments and 
simulations respectively. 

We investigated the nonlinear propagation for 2 < 
\S\ /r < 12 and 25 ^iW < P < 1 mW. A general fea- 
ture is that the nonlinear effects, heralded by diverging 
behaviors for opposite Ss, start to develop roughly 2 mm 
before the peak density is reached (see Fig. 2). The sub- 
sequent behavior for z > can then be characterized by 
the "divergence" a = d(Ay)/dz. In Fig. 3 we report a 
vs the laser power, for |(5| /Y = 4. Circles and squares 
correspond to S > and < 0, respectively (filled and 
empty symbols for experiment, respectively, simulation). 
The horizontal solid line represents the divergence of the 
incident Gaussian beam. The dotted line a = corre- 
sponds to self-trapping. For a < the beam is strongly 
focused inside the cloud, while for a > the beam di- 
verges, either due to insufficient self-focusing (6 > 0) or 
self-defocusing (S < 0). The comparison between exper- 
iment and numerics is quite satisfying. In Fig. 4a) and 
b) we show numerical plots of the propagating light in- 
tensity in the y — z plane for the parameters of Fig. 2. 
The contour lines represent the density of the cold atom 
cloud. In the self- focusing case b), a stable self-trapping 
of the beam is observed in the region shown by the gray 
arrow. 

In conclusion, we have provided an experimental evi- 
dence of light self-trapping in a large cloud of cold atoms. 
Numerical simulations based on the microscopic nonlin- 
ear atomic susceptibility are in good agreement with the 
experiment. The observation of such a solitonic behavior 
paves the way for the study of nonlinear physics with 
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Fig. 4. Numerical plots of the light intensity in the y — 
z plane for the parameters of Fig. 2: a) self-defocusing 
and b) self-focusing case; the contour lines represent the 
density of the cold atom cloud in cm~ 3 . The gray arrow 
indicates the region where the beam waist is constant. 



light in cold atoms, a medium with unique possibilities 
in terms of tunability. In particular, this work consti- 
tutes a first step towards the study of optical localized 
structures in cold atoms. In contrast with the situation 
in hot vapors [11], one expects here a strong interaction 
between the localized structures and the atomic external 
degrees of freedom, with a consequent complex dynam- 
ics [12]. 
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